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Abstract. We review the actual state in the description of the NN interaction 
by means of chiral constituent quark models. We present a series of relevant 
features that are nicely explained within the quark model framework. 

1 Introduction 

At the end of the 70 's the potentialities of quark model calculations in low- 
energy hadron physics were established. In a series of pioneering works, Isgur 
and Karl Q performed a quite successful exhaustive study of the baryon spec- 
tra and at the same time the short-range repulsion of the nuclear force was 
qualitatively explained Q . These ideas encouraged several groups to undertake 
the ambitious project of trying to understand the NN interaction in terms of 
quark degrees of freedom. 

Several authors ||, ||, |j| found a quantitative explanation of the NN short- 
range repulsion based on dynamical effects induced by the one-gluon exchange 
(OGE). One would say that there exits no repulsive potential in the naive 
sense. The color-magnetic OGE interaction likes to have as many symmetric 
pairs in the color-spin space as possible. Due to the connection of the spin- 
color space with the orbital space as a consequence of the Pauli principle the 
orbital symmetry [42] o is preferred. But this symmetry has to have a zero in 
the S-wave. This zero yields for the phase shift the same results as a hard or a 
soft core. 

Based on this success simple quark models of the nuclear force were con- 
structed. The so-called hybrid quark models Q contained three main ingredi- 
ents: (i) A short range part given by the quark-exchange interaction, (ii) a long 
range part by the one-pion exchange (OPE) potential, and (iii) a medium range 
part by a phenomenological potential or a two-pion exchange potential. These 
meson exchange potentials were considered to simulate the meson cloud sur- 



2 



rounding the quark core. These models allowed to reproduce the NN scattering 
and bound state data. The next conceptual step in order to understand the NN 
interaction only based on quark degrees of freedom was the consideration of 
(qq) and (qq) 2 excitations generated by off-shell terms of the Fermi-Breit quark- 
gluon interaction Although the obtained potential acquired an attractive 
part in the 0.8—1.5 fen range, it was too weak to bind the deuteron or to fit 
the extreme low-energy S-wave scattering. 

In the middle 80's quark antisymmetry effects on the OPE started to be 
analyzed Jl}. A similar NN short-range repulsion to the one provided by the 
OGE was obtained. A first quark-model calculation based on gluon and pion 
exchange at the level of quarks was performed ||. As the OPE contributes 
to the A — N mass difference, the rather big quark-gluon coupling constant 
used in the previous models (a s > 1) got nicely reduced to values around 
0.5. However, it was not obtained enough intermediate range attraction. To 
avoid this problem, a phenomenological scalar potential at baryonic level was 
introduced, returning on this way to a hybrid model. Besides, its coupling 
constant was fitted to a different value for NN S'-waves, than for higher angular 
momentum partial waves or for the deuteron 0] . 



2 The Chiral Constituent Quark Model 

In the early 90 's the constituent quark mass was related to the breaking of 
chiral symmetry [|[ [ic[| . The underlying idea is that the constituent quark 
mass is a consequence of the spontaneous chiral symmetry breaking. Then, be- 
tween the chiral symmetry breaking scale (vIcsb ~ 1 GeV) and the confinement 
scale, (Aq, ~ 0.2 GeV) QCD may be simulated in terms of an effective theory 
of constituent quarks interacting through Goldstone modes. Specifically, in a 
Nambu-Goldstonc realization of chiral symmetry within the linear sigma model, 
there appear two Goldstone boson fields: the pion and the sigma. The chiral 
constituent quark model incorporates OPE and one-sigma exchange (OSE) in- 
teractions of the form, 
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A being a cutoff parameter and Y(x) and H(x) the standard Yukawa functions. 
Perturbative features of QCD are incorporated through the OGE potential, 
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where the A's stand for the color SU(3) matrices. Finally, a phenomenological 
confining potential is introduced. Its detailed radial structure, being funda- 
mental to study the hadron spectra, is expected to play a minor role for the 
two-baryon interaction. To be consistent with meson spectroscopy it can be 
taken to be linear 

Vcon(7«) = ~a c Ai • \j nj , (4) 

though for simplicity harmonic oscillator quark wave functions will be used for 
the NN interaction. 

This model is able to reproduce quite nicely the NN scattering and bound 
state data Jl^, |ll[ . Typical values for the parameters are given in Table [l]. Other 
quark-model approaches to the NN interaction are found in the literature [lL2| . 



Table 1. Quark-model parameters. 
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3 Several Reasons to Use Quark Models 

3.1 Medium Range Attraction 

A new feature of the chiral constituent quark model is that through the OSE 
potential incorporates in a natural way the NN medium range attraction with- 
out any additional parameter. This fact is illustrated in Fig. |l| where we com- 
pare the scalar potential obtained at quark level (solid line) with the baryonic 
parametrization used in ref. @] (dashed line) to fit the experimental data. While 
at the baryonic level two different coupling constants are used: ff^ NN /4rK = 3.7 
for (S,T)=(0,1) and & NN /4>ir = 1.9 for (S,T)=(1,1), the quark model result is 
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Figure 1. iVTV a exchange potential 
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obtained in both cases with the same parameters fixed from the irqq coupling 
constant. In Fig. || we compare the phase shifts for both cases and we can see 
how the predictions are similar. 

3.2 NN Short-Range Repulsion 

The NN short-range behavior has been explained in terms of an interplay 
between quark antisymmetry and a particular quark dynamics. In a group 
theory language, a completely antisymmetric six quark state asymptotically 
describing two free nucleons in a relative S-wave is given, 
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the subindex C, O or ST indicating the color, orbital or spin-isospin representa- 
tion respectively. When the distance between the two clusters goes to zero the 
{42} configuration increases in energy. According to our harmonic oscillator 
wave function ansatz for the quarks, the excitation energy is given by, 



A ho = A [£ { 42} -£{6}] ho = 2huj = 



2h 2 
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Let us now consider other interactions. Assume first that quarks interact via 
the OGE used in ref. § (a s = 1.39, m q = 300 MeV and b = 0.6 fm). We 
can estimate the energy contribution of the OGE calculating the matrix el- 
ement, < {} | J2i<j Voge(t ij) | {} > where {} stands for the <^ 42 } and 
^{gj configurations of Eq. (g). As we are interested in the NN interaction 
energy we should subtract twice the corresponding nucleon self-energy given 
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by < N | J2i<j VoGE(nj) I N >. For (S,T)=(1,0), the matrix element of the 
color magnetic potential for the {6} configuration is 48.83 MeV, while for the 
{42} configuration is — 235.63 MeV. Since in the energy difference between 
the {42} and {6} configurations the subtracted self-energies cancel, we have 
zAoge = A [£'{42} — -^{6}] ge = — 284.46 MeV. This difference compensates 
the harmonic oscillator energy difference, Ziho = 240.36 MeV, the [42] o orbital 
symmetry becoming the lowest in energy. As a consequence of this configu- 
ration mixing the S-wave NN relative motion wave function shows an energy 
independent node, which translates into a hard-core behavior in the phase shift. 

In the case of the chiral constituent quark model the value of a s , which 
drives the OGE energy gap between the {42} and {6} configurations, is sig- 
nificantly reduced due to the pseudoscalar contribution to the A — N mass 
difference and correspondingly the mixing effect. If we recalculate the contri- 
bution of the OGE with the parameters of Table [l], we observe again a de- 
crease of the energy difference between the {42} and the {6} configurations 
^OGE = —140.99 MeV, but it is much smaller than the harmonic oscillator 
energy difference Z\ho = 309.10 MeV and it is not sufficient to give rise to a 
significant configuration mixing. It is then convenient to deep in the under- 
standing of the repulsion mechanism. 

A similar kind of analysis can be carried out for the OPE and OSE po- 
tentials. For (S,T)=(1,0), differently than in the OGE case, the OPE has the 
same sign for both configurations. Regarding the energy difference one gets 
^ope = -68.62 MeV. Concerning the OSE one has Z\ se = 133.05 MeV. 
Therefore, in the chiral constituent quark model three effects conspire against 
the symmetry mixing. First, the small value of a s . Second, the cancellation 
between the OPE contributions to the {42} and {6} configurations, and finally 
the contrary effect of the OSE potential. Putting all the contributions together 
one obtains, A [£{42} - ^{6}] oge+ope+ose = -76.56 MeV, which is much 
smaller than the harmonic oscillator energy difference Z\ho = 309.10 MeV. 
Then, one must conclude that in the chiral constituent quark model there is 
not enough configuration mixing to account for the hard-core of the NN inter- 
action as a node produced by the [42] o orbital symmetry. 

To look for the origin of this repulsive character of the interaction one 
should go beyond the energy difference between the symmetries and calcu- 
late the specific contribution of the interaction for each symmetry. In or- 
der to see the repulsive or attractive character of every term of the poten- 
tial in both configurations one has to subtract twice the corresponding nu- 
cleon self-energy, given by £oge = —72.63 MeV, £ope = —311.40 MeV, and 
£ OSE = -66.90 MeV. One then obtains, e°° e+ope+ose = 455.72 MeV and 

E°^ E+OPE+OSE = 379.16 MeV. It is then clear that the chiral potential pro- 
duces strong repulsion in both configurations, {42} and {6}. The OPE hardly 
contributes to the mixing between them, however in both symmetries it pro- 
duces strong repulsion. The OGE reduces the energy difference whereas the 
OSE increases this energy difference and produces an additional attraction in 
both configurations. The whole effect is an energy difference of about the same 
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value obtained with the harmonic oscillator model but with an additional re- 
pulsion in both symmetries originated mainly by quark antisymmetry on the 
OPE. Therefore, in this type of models the NN S-wave behavior should be 
attributed to the strong repulsion generated by the OPE in the [6]o orbital 
configuration. 



3.3 Universality 

A main feature of the quark treatment is its universality in the sense that all 
the baryon-baryon interactions are treated on an equal footing. Moreover, once 
the model parameters are fixed from NN data there are no free parameters 
for any other case. This allows a microscopic understanding and connection 
of the different baryon-baryon interactions that is beyond the scope of any 
analysis based only on effective hadronic degrees of freedom. We will illus- 
trate our discussion by means of the recently calculated NN — * TV TV* (1440) 
transition potential In particular we will determine the 7riViV*(1440) and 
crNN* (1440) coupling constants. 

For this purpose, let us realize that asymptotically (R > 4 fm) the OSE and 
OPE potentials have at the baryon level the same spin-isospin structure than at 
the quark level. Hence we can parametrize the asymptotic central interactions 
as, 
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where gi stands for the coupling constants at the baryon level and M r is the 
reduced mass of the iViV*(1440) system. 

By comparing these baryonic potentials with the asymptotic behavior of 
the OPE and OSE obtained from the quark calculation we can extract the 
7riViV*(1440) and <riViV*(1440) coupling constants in terms of the elementary 
nqq coupling constant and the one-baryon model dependent structure. The sign 
obtained for the meson- TV TV* (1440) coupling constants comes determined by 
the arbitrarily chosen relative sign between the N and TV* (1440) wave functions. 
Only the ratios between the ?riViV*(1440) and aNN* (1440) would be free of 
this uncertainty. 

To get g (1440) / 'y/far we turn to our numerical results for the 1 So OPE 
potential, and fit its asymptotic behavior to Eq. (||). We obtain 

toNNgyNwrn a 2 3 (10) 
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i.e. g-KNN*(iAAo)l \P^k = —0.94. As explained above only the absolute value of 
this coupling constant is well defined. The coupling scheme dependence can be 
explicitly eliminated if we compare 3 7r A'Ar*(i44o) with g^NN extracted from the 
NN — > NN potential within the same quark model approximation. We get 



|<?7riVAr*(1440)/#7rArjv| = 0.25. 

By proceeding in the same way for the scalar potential we can write 

\gaNN*(144Q)/gaNN \ = 0.47. 



(11) 



(12) 



The ratio given in Eq. (jll|) is similar to that obtained in ref. ]14| . Nonetheless 
one can find values for f n n jv* (1440) ranging between 0.27—0.47 coming from 
different experimental analyses. Regarding the ratio obtained in Eq. (12), our 



result agrees quite well with the only experimental available result, obtained 
in ref. [[l5| from the fit of the cross section of the isoscalar Roper excitation in 
p(a, a') in the 10—15 GeV region, where a value of 0.48 is given. Finally, we 
can give a very definitive prediction of the magnitude and sign of the ratio of 
the two ratios, 



57rWA r *(1440)/ 'QttNN — 0.53 g a N N * (1440) /g a NN ■ 



(13) 



3.4 Pauli Blocking 

The Pauli principle plays an essential role to explain the atomic and nuclear 
structure. For example, the strong repulsion existing in the ctu system origi- 
nates from the nucleonic substructure of the a particle and not from a particular 
dynamics. The hard-core appears because the saturation of nuclear spin-isospin 
degrees of freedom allows only four nuclcons in the orbital ground state. We 
shall refer to this saturation phenomenon indicated by the existence of forbid- 
den states at the underlying structure level as Pauli blocking. At the level of 
subnuclear physics the elementary constituents, the quarks, are fermions so a 
quark Pauli principle is active. 

For identical baryons the symmetrization postulate at the baryonic level 
gives rise to selection rules on the state of the compound system. At the quark 
level, assuming SU(2) flavor symmetry, this selection rule appears as a partic- 
ular case of the action of the quark antisymmetrizer, A, on a two cluster state 
from six identical quarks. Taking into account that any two-baryon state (LST) 
can be decomposed in a symmetric plus an antisymmetric part under the ex- 
change of the baryon quantum numbers, one can write for a definite symmetry 
(specified by /) 
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where B\ and B 2 denote the baryons, and — 4p, y are tne positions of the 
clusters. The antisymmetrizer can be decomposed as, 

A=±(l-9P aa )(l-P), (15) 

where P 3 6 exchanges quarks 3 and 6, and V exchanges the two clusters. The 1 — 
V factor implies L + Si + S 2 - S + Ti + T 2 - T+/ = odd. For non-identical 
baryons this relation indicates the symmetry associated to a given set of values 
of (LST). For identical baryons, B± = B2, / = even (in order to have a no 
vanishing wave function) and L + S + T = odd, reproducing the well-known 
selection rule. Certainly the effect of the quark substructure goes beyond the 
(1 — V) factor in the antisymmetrizer, being also included through the terms 
containing P 36 . Let us study the effects of these quark exchanges in a system of 
two non-identical baryons, in particular the NA system. Assuming the three- 
quark cluster wave function at position S given by 

3 

<P B (123; S) = H - S) ® [3] ST ® [l]c , (16) 

i=l 

where rj 0s {x) is a gaussian with parameter 6, and [3]st ([1]c) stands for the 
spin-isospin (color) representation, the norm of the projected ^jfj state reads, 

NjfJ(R) - N?(R) - C(S, T, /) Nr(R) , (17) 

where N'^(R) [Nf^(R)] refers to the direct [exchange] radial contribution com- 
ing from the l[P3e] term in the antisymmetrizer. The spin-isospin coefficient 
C(S,T,f) is given by, 

C(S,T,f) = 3 | ST (iV(123)Z\(456) | P 3 S 6 T | 7V(123)Z\(456)) ST + 

(-1)' st(^V(123)Z\(456) | P 3 5 6 t | Z\(123)7V(456)) S t} .(18) 
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Its value for the N A S'-waves is given in Table I (the NN case is also shown for 
comparison). In the R — ► limit and to the dominant order Eq. (fL^) transforms 



Table 2. C(S, T, /) coefficients for the NN and NA S-waves. 
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Then, for 3 L = C(S,T,f) the overlapping of the two-cluster wave function 
behaves as R 2L+A instead of the centrifugal barrier behavior R 2L , indicating 
that quark Pauli blocking occurs (in other words, a node appears in the rel- 
ative wave function). This turns out to be the case (see Table ||) for the NA 
partial waves 3 Si(T = 1) and 5 S , 2(T = 2) giving rise to forbidden states (the 
[6]o orbital configuration for L=0). If we now calculate the NA potential, we 
observe (solid line) the important difference between the potential in a Pauli 
blocked and a non Pauli blocked channel (see Fig. ||). At the same time we 
have represented by the dashed line a NA potential obtained from a direct 
scaling of the NN interaction jl6| and we see how the difference between the 
two channels is just a spin-isospin factor. As shown in Fig. ||, a separable NA 
potential reproducing the phase shifts obtained from the experimental ird elas- 
tic cross section presents the strong repulsion predicted by the quark model in 
the 3 S , i(T = 1) partial wave @- 

A similar analysis for the NN system does not show any forbidden con- 
figuration, see Table |[ As we have seen, the NN short-range repulsion in the 
S- waves involves dynamical effects. 

3.5 Connection to Other Systems: The Baryon Spectrum 

The chiral constituent quark model allows also to study the baryon spec- 
trum in a completely parameter-free way. We have calculated the nonstrange 
baryon spectrum within a Faddeev approach. We included all the configura- 
tions (£, A, s, t) (£ is the orbital angular momentum of a pair, A is the orbital 
angular momentum between the pair and the third particle, while s and t are 
the spin and isospin of the pair) with t and A up to 5. 

For this purpose, the delta function of the OGE interaction has to be reg- 
ularized. By choosing an exponential regularization 



*(»■<. 



1 



-Uj/ro 



Airr 2 , 



(20) 



with ro = 0.25 fm (a typical value for spectroscopic models) we obtain the N 
and A spectra shown in Fig. The predicted spectrum is reasonable except 
for the fact that all the negative parity states are around 100 MeV below the 
experimental data. The relative position of the N(1AA0) and the A(1535) is not 
very much affected by the modification of a s (see Fig. 3 of ref. |l8) ) , indicating 
that the correct level ordering is not connected with the strength of OGE. 

A possible mechanism responsible for the inversion of these states can be 
obtained through the pseudoscalar interaction. To illustrate this point we have 
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Figure 5. Relative energy nucleon and A spectra 

calculated the energy of the iV(1440) and AT(1535) increasing the contribution 
of the pseudoscalar interaction by letting the cutoff parameter of the OPE 
to increase. The results are shown in Fig. |(| As can be seen, the inversion of 
the ordering between the positive and negative parity states can be achieved if 
A n becomes sufficiently large (around 7 fm _1 for the set of parameters of Table 
. Such a model would be incompatible with the understanding of the basic 
features of the two- nucleon system [fill . 

4 Summary 

In this work we have emphasized several physical aspects of the NN inter- 
action that support the results obtained within the chiral constituent quark 
model. In particular, medium range attraction is satisfactorily explained and 
the NN short-range repulsion is understood through antisymmetry effects on 
the OPE interaction. Actually, antisymmetry effects at the level of Goldstone 
bosons appear as a basic tool to understand the NN phase shifts. Extension of 
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Figure 6. iV(1440) and iV(1535) masses as a function of A w 
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the model to the general baryon-baryon case is straightforward and parameter- 
free. From it, non-strange meson-baryon-baryon coupling constants can be cal- 
culated. Pauli blocking effects manifest in different partial waves of several 
systems and find support in indirect experimental data. Such strong repulsion 
cannot be understood in a meson-exchange model. 

Finally, the chiral constituent quark model is able to generate a quite rea- 
sonable description of the baryon spectrum with a set of parameters that allows 
to understand the NN phenomenology. Therefore, the validity of a model where 
the OGE is combined with Goldstone-boson exchanges is out of question, and 
it presents advantages with respect to models based only on Goldstone-boson 
exchanges, as has been recently emphasized in several works ||l9| , |20| . 
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